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摘要 : 以 河北 省 中 南部 面积 为 93 539 平方 公里 ( 约 占 河北 总 面积 的 51.5 和 % ) 的 区 域 作 为 研究 区 ,应 用 地 理 信 息 系 统 技术 和 地 统 
计 学 方法 对 1980 ~ 1997 年 的 来 自 36 个 监测 点 的 二 代 棉 铃 虫 卵 的 密度 数据 进行 了 空间 结构 分 析 和 空间 分 布 模拟 。 结 果 表 明 旷 
在 所 研究 的 尺度 上 存在 空间 相关 性 ,但 各 年 空间 相关 的 程度 均 较 弱 且 随 年 份 和 地 域 的 不 同 而 不 同 。 空 间 变 异 成 分 所 占 比例 
的 变化 范围 为 3.54% ~ 61.67% ,其 相关 程 从 115.92 km 到 487.08 km。 时 间 序 列 分 析 表 明 , 卵 的 发 生 数 量 始 终 处 于 波动 状态 ， 
其 中 1980 ~ 1991 年 维持 在 低 水 平 的 波动 状态 ,其 发 生 的 平均 密度 为 330 粒 / 白 株 , 变化 范围 为 122~ 706 粒 / 百 株 ; 而 1992 ~ 1997 
年 在 高 水 平 波 动 , 其 平均 密度 为 1 877 粒 / 百 株 ,变化 范围 为 264 ~4 023 粒 / 百 株 。 空 间 模 拟 结 果 表 明 , 各 年 的 卵 在 空间 均 呈 现 
聚集 分 布 的 格局 ,但 聚集 程度 依据 其 空间 相关 成 分 所 占 比 例 的 不 同 而 不 同 。 从 发 生 的 地 域 特 点 看 ,从 东北 至 西南 , 卵 的 发 生 
程度 逐渐 趋 重 , 其 中 西部 的 阜 平 县 .曲阳 县 、 行 唐 县 、 有 灵 寿 县 .平山 县 等 地 以 及 南部 的 基 郸 市 、 磁 县 、 临 漳 县 、 魏 县 等 地 发 生 最 重 ， 
而 东北 的 遵 化 市 、 玉 田 县 .丰润 县 .唐山 市 、 丰 南 县 等 地 最 轻 。 卵 的 密度 的 各 个 等 级 不 仅 在 地 域 上 存在 差异 ,而且 随 时 间 发 生 
变化 ,从 整体 发 生 情 况 看 ,1992 年 发 生 最 重 , 其 次 为 1993 和 1995 年 ,1981 和 1990 年 发 生 较 轻 , 1987 年 发 生 最 轻 , 1980 年 代 的 发 
生 数 量 显著 低 于 1990 年 代 , 这 与 时 间 序 列 的 分 析 结 果 一 致 。 棉 铃 虫 1990 年 代 的 发 生 偏重 主要 是 气温 持续 偏 高 .冬季 变 暖 所 
致 。 而 研究 区 棉铃 虫 重 发 生地 位 于 西部 及 西南 部 分 地 区 ,该 地 的 主要 特点 是 多 山 、 植 被 丰富 .小麦 种 植 比例 高 .6 月 降雨 量 偏 
少 , 根 据 相 关 研究 ,多 样 性 的 寄主 和 干旱 的 气候 条 件 有 利于 棉铃 虫 的 大 发 生 。 

关键 词 : 棉铃 虫 ; 卵 ; 地 理 信息 系统 ;地 统计 学 ;空间 分 布 ;动态 

中 图 分 类 号 : 0968 ”文献 标识 码 : A ”文章 编号 : 0454-6296(2004)01-0033-08 

Analysis on spatial distribution and dynamics of Helicoverpa armigera \ Hiibner ) 


eggs, based on GIS and GS 

WANG Zheng-Jun; LI Dian-Mo”* , XIE Bao-Yu (State Key Laboratory of Integrated Management of Pest Insects and Ro- 
dents，Institute of Zoology，Chinese Academy of Sciences, Beijing 100080, China) 

Abstract: In this study, Geographic Information System《GIS) technology and geostatistics (CGS) methods were used to 
analyze the spatial structure and simulate the spatial distribution with the historical data of the 2nd generation eggs densi- 
ties of Helicoverpa armigera (Hiibner)}, which were collected from 36 monitoring points in Hebei Province during 1980 to 
1997. The study areas lie in central and southerm part of Hebei Province, and amounted to 93 539 km ，about 51.5% of 
the whole provincial area. Results showed that the eggs were correlated spatially in the study areas, but degrees of corre- 
lation were all lower and different from year to year and area to area. Proportion of spatial variation changed from 3.54% 
to 61.67%; the distance range of spatial dependence varied from 115 .92 to 487.08 km. The time series analysis showed 
that the quantity of H. armigera was always in fluctuation during 1980 to 1997. During the period of 1980 to 1991, the 
quantity was low with average density of 330 eggs per a hundred plants; and its range was 122 - 706 eggs per a hundred 
plants. However, the egg quantity became high remarkably during 1992 to 1997 with average density of ] 877 eggs per a 
hundred plants, and its range was 264 - 4 023 eggs per a hundred plants Based on the outbreak map of the 2nd genera- 
tion of H. armigera, it was found 也 at the spatial distribution patterns were aggregated in all years, but the collective de- 
gree was changed with time and the proportion of the spatial correlation component. The density grades of egg distribution 


were high in the southwest, but low in the northeast in the study areas. From the northeast to the southwest, the densities 
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of eggs increased gradually. The highest grades of density were found in the areas of Fuping’; Quyang’ Xingtang, Ling- 


shou, Pingshan, Handan; Cixian; Linzhang and Weixian; and the lowest in Zunhua, Yutian, Fengrun, Tangshan and 


Fengnan. The density of eggs was not only related to location, but also to time. On the whole，outbreak was the heaviest 
in 1992 and heavier in 1993 and 1995, while moderate in the 1981 and 1990 and the lightest ln 1987. The outbreak 


number in the 1980s was significant lower 也 an that in the 1990s. The results were consistent with that of the time series 


analysis. Based on the pertinent Tepotts and our study, it was suggested that the higher air temperature, especially warm- 


er winter was the main factor which triggered outbreak of the H. armigera ln the 1990s. It was also suggested that in 


west and southwest outbreak areas, numerous mountains, abundant vegetation, higher wheat cultivation proportion, and 


lesser rainfall in June, should be important factors leading to the outbreak of H. armigera. 
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了 解 害虫 的 区 域 性 空间 分 布 格局 及 动态 是 进行 
害虫 宏观 管理 的 前 提 和 基础 。 只 有 掌握 害虫 的 空间 
动态 和 时 间 动 态 , 才能 进一步 地 揭示 其 暴发 的 机 制 
并 对 其 未 来 的 发 生 进行 预测 。 

目前 ,对 害虫 区 域 性 空间 分 布 及 动态 的 研究 仍 
然 是 一 个 薄弱 环节 。 这 表现 在 两 个 方面 : 首先 , 许 
多 研究 仍 局 限于 时 间 序 列 而 忽视 了 空间 因子 对 种 群 
动态 的 影响 ( 程 新 意 和 杨 崇 瑞 , 1992; 黄 荣 花 等 ， 
1995; 宫 锡 鸿 等 , 1997)。 例 如 各 地 植保 站 在 进行 监 
测 活动 时 普遍 存在 抽样 点 偏 少 , 以 点 盖 面 现象 ,而 昆 
忠生 存 的 环境 往往 是 复杂 的 、 异 质 性 的 ,用 局 域 的 发 
生来 代表 整体 的 情况 其 结果 必然 带 来 监测 和 预测 上 
的 偏差 。 此 外 ,各 地 站 点 间 在 害虫 监测 上 还 缺少 横 
回 的 联系 和 沟通 , 一般 多 关心 本 地 害虫 的 活动 而 急 
略 了 异地 害虫 的 发 生 对 当地 的 影响 。 其 次 ,就 区 域 
性 害虫 空间 动态 而 言 , 研究 比较 多 的 是 林业 害虫 如 
舞 毒 蛾 ( Hohn et ail.,，1993; Zhou et al.,， 1995; Lieb- 
hold et al.，1998), 水稻 害 忠 如 福 飞 恒 ( 王 海 扣 等 ， 
1998) 等 。 由 于 害虫 数据 的 限制 及 缺乏 有 效 的 分 析 
工具 ,国内 在 这 方面 的 研究 仍然 很 不 成 熟 , 仅 涉 及 很 
少 的 害虫 如 稻 飞 剧 、 二 化 蜡 等 ( 王 正 军 等 ,2000; 周强 
等 ,2001,2003), 一 些 涉及 覃 铃 虫 的 研究 也 仅见 于 小 
尺度 (如 田 块 水 平 ) 和 空间 结构 分 析 ( 葛 绍 奎 等 ， 
2000; 陆 永 跃 和 深 广 文 ,2002), 对 大 尺度 的 棉铃 虫 空 
间 分 布 和 动态 的 研究 还 未 见报 道 。 因 此 如 何 利 用 有 
限 的 历史 数据 .合理 的 分 析 方 法 来 研究 和 分 析 害 上 虫 
区 域 性 暴发 成 灾 的 规律 性 是 目前 我 国 昆 虫 生 态 学 工 
作者 需要 关注 的 课题 。 

棉铃 虫 是 我 国 主要 棉花 害虫 之 一 ( 翟 连 荣 等 ， 
1992; 吴 坤 君 和 获 佩 瑜 ,2001), 而 河北 省 是 我 国 重要 
的 覃 花 产 区 ,棉铃 虫 在 该 省 的 发 生 世 代为 三 代 或 四 
代 , 其 中 以 二 代为 害 最 重 , 为 主 害 代 。 二 代 棉 铃 虫 在 
华北 棉 区 主要 集中 在 棉田 内 发 生 为 害 , 而 三 代 棉 铃 


虫 则 分 散在 棉田 、 玉 米田 、 花 生 田 内 (6 丁 岩 钦 和 张 占 
川 ,1994>。 由 于 二 代 卵 粒 基 数 的 大 小 直接 影响 全 年 
的 发 生 情况 ,加 之 卵 在 空间 分 布 上 的 相对 稳定 性 ， 
此 ,我 们 以 河北 省 棉铃 虫 二 代 卵 作为 研究 对 象 ,对 其 
在 历史 上 的 发 生 情况 进行 系统 分 析 , 其 主要 目的 有 
三 个 ; 1 明确 覃 铃 虫 卵 在 区 域 性 尺度 是 否 存在 空间 
相关 关系 ;2) 通 过 地 统计 学 的 方法 获得 各 年 卵 的 衬 
间 分 布地 图 ;37 明确 二 代 棉 铃 虫 卵 在 空间 的 分 布 格 
局 及 其 动态 规律 。 


1 资料 和 方法 


1.1 研究 区 域 与 数据 

以 河北 省 中 南部 地 区 (面积 为 93 539 平方 公 
里 , 约 占 河北 总 面积 的 $1.5% ) 作 为 研究 区 域 , 共 涉 
及 的 抽样 点 为 36 个 (图 1)。 选 择 这 个 区 域 主要 基 
于 以 下 几 方 面 的 考虑 ; 1) 该 地 为 河北 省 棉花 种 植 的 
主要 区 域 ,北部 高 纬度 地 区 棉花 种 植 比 例 很 小 , 明确 
这 一 地 区 的 棉铃 忠 种 群 动态 规律 可 以 反映 全 省 的 情 
况 ;2) 由 于 1980 ~ 1997 年 这 18 年 的 抽样 点 大 都 位 于 
这 一 地 区 ,对 这 些 抽 样 点 进行 空间 模拟 可 以 产生 比 
较 合理 的 插值 结果 。 

研究 使 用 的 数据 包括 : 研究 区 边界 图 (1: 
1 000 000); 研 究 区 数字 高 程 模型 DEM(1:250 000)。 
各 抽样 点 5 月 和 6 月 各 旬 的 平均 气温 、 旬 降雨 量 ;5 
月 和 6 月 的 月 平均 气温 、 月 降雨 总 量 、 温 雨 系数 、 温 
湿 系 数 、 降 雨 强度 以 及 高 程 数据 、 各 年 小 麦 和 棉花 的 
种 植 面积 比例 等 。 
1.2 研究 方法 
1.2.1 空间 结构 分 析 和 空间 分 布 模拟 : 区 域 性 变 
量 在 空间 的 分 布 主要 由 两 部 分 组 成 , 即 空间 结构 性 
成 分 (用 偏 基 台 值 与 基 台 值 的 比 表 示 其 所 占 比 例 的 
大 小 ) 和 随机 性 成 分 (用 块 金 系 数 与 基 台 值 之 比 表示 
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图 1 研究 区 和 抽样 点 所 在 位 置 
Fig.1 The locations of the study areas and sampling points 


其 所 占 比 例 的 大 小 )。 其 中 空间 结构 性 成 分 由 空间 
相关 而 产生 的 ,而 随机 性 成 分 来 源 于 随机 性 误差 。 
空间 结构 分 析 应 用 协 方差 或 半 方 差 销 数 进行 ,对 于 
卵 的 空间 结构 ,本 研究 使 用 半 方 差 图 和 协 方差 图 两 
种 方法 进行 模拟 ,模拟 模型 主要 有 三 种 ,包括 球 型 模 
型 指数 模型 和 高 斯 模型 ,具体 选择 时 主要 根据 模拟 
误差 最 小 的 原则 。 

空间 模拟 所 用 插值 方法 包括 普通 克 立 格 法 (Or- 
dinary kriging) ,普通 协同 殉 立 格 法 (Ordinary Cokrig- 
ing) 和 简单 克 立 格 法 (Simple kriging) 三 种 。 参与 协 
方差 分 析 和 协同 克 立 格 插值 的 变量 有 :5 月 和 6 月 
各 名 的 平均 气温 、 旬 降雨 量 ;$ 月 和 6 月 的 月 平均 气 
温 ,月 降雨 总 量 . 温 雨 系数 , 温 湿 系数 .降雨 强度 以 及 
高 程 数据 .各 年 小 麦 和 棉花 的 种 植 面积 比例 等 。 空 
间 模 拟 的 结果 将 用 于 空间 估 值 ,空间 估 值 采用 邻 域 
搜索 法 。 

各 年 卵 分 布 的 模拟 地 图 的 生成 方法 按 以 下 方式 
确定 : 首先 用 Ordinary kriging, 如 果 结 果 误 差 足 够 
小 , 即 接 近 于 零 ,将 不 再 用 Ordinary Cokriging 插值 ; 
如 果 误 差 较 大 , 则 与 Ordinary Cokriging 方法 进行 比 
较 , 选 择 误差 较 小 者 作为 最 终 的 插值 方法 。 应 用 
Ordinary Cokriging 方法 时 ,以 二 代 卵 的 密度 为 主 变 
量 ,其 它 因子 作为 协 变量 进行 插值 ,通过 对 各 个 因子 
逐一 进行 插值 比较 ,选择 插值 误差 最 小 的 那个 因子 
作为 协 变量 生成 最 终 的 卵 分 布 图 -。 

为 了 定量 分 析 卵 的 发 生 , 将 各 年 的 插值 图 分 为 
5 个 等 级 , 即 轻 发 生 ( < 50 粒 / 百 株 )、 中 偏 轻 (51 ~ 
150 粒 / 百 株 ) .中 度 发 生 (151 ~ 500 粒 / 百 株 ) .中 偏重 
(501 ~ 1 500 粒 / 百 株 ) 和 大 发 生 ( > 1 500 粒 / 白 株 )。 
根据 各 年 的 发 生 分 布 图 ,统计 各 发 生 等 级 的 面积 和 
比例 。 

空间 结构 分 析 和 分 布 模拟 是 在 地 理 信息 系统 软 
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件 ArcGIS 8.1(ESRI Inc. 2001) 和 ArcView GIS 3.2 
(ESRI Inc. 1999) 的 支持 下 进行 的 。 
1.2.2 交叉 验证 : 所 谓 交 义 验 证 是 指 移 去 一 个 数 
据点 ,用 其 它 数 据 来 估计 该 点 以 检验 估 值 精度 的 方 
法 。 对 于 一 个 较 好 的 模型 ,应 该 满足 : 1) 其 平均 预 
测 误差 应 该 接近 于 0;2) 均 方 根 误差 应 尽 可 能 的 小 ; 
3) 平 均 标准 差 应 该 接近 于 0;4) 平 均 克 立 格 标准 差 
应 尽 可 能 的 小 ;5) 均 方 根 标准 差 应 该 接近 于 1;6) 均 
方 根 误差 与 平均 克 立 格 标准 差 之 间 差 的 绝对 值 也 要 
尽 可 能 的 小 。 模 拟 结 果 将 根据 以 上 标准 来 比较 、 判 
断 和 选择 。 

选择 最 优 模型 的 路 线 为 : 1) 首 先 选择 满足 平均 
预测 误差 尽 可 能 的 接近 于 0;2) 选 择 平均 标准 差 尽 
可 能 的 接近 于 0;3) 从 中 选择 满足 均 方 根 标准 差 接 
近 于 1;4) 最 后 再 从 中 选择 出 均 方 根 误 差 与 平均 克 
立 格 标准 差 之 间 差 的 绝对 值 最 小 的 方法 ,也 就 是 最 
优 方法 。 


2 结果 与 分 析 


2.1 卵 的 空间 结构 分 析 

各 年 的 插值 方法 和 相应 的 模型 参数 如 表 1 所 
示 , 对 于 普通 协同 克 立 格 方法 ,我 们 仅 选 择 二 代 卵 的 
密度 的 空间 结构 参数 ,不 包括 卵 与 其 它 协 变量 及 协 
变量 本 身 的 变异 情况 ,空间 变异 成 分 和 随机 程度 也 
仅 针对 卵 本 身 而 言 。 

根据 表 1 可 以 看 出 : 1) 二 代 棉 铃 虫 卵 在 各 年 都 
存在 一 定 的 空间 相关 性 ,空间 相关 距离 (相关 程 ) 为 
115.92 ~ 487.08 km;2) 卵 的 密度 的 空间 变异 成 分 在 
各 年 都 比较 低 , 变 化 范围 为 3.54% ~61.67%-. 
2.2 各 年 卵 的 空间 分 布 图 的 生成 

应 用 地 统计 学 的 插值 方法 和 模型 ,在 空间 相关 
性 分 析 的 基础 上 ,对 河北 省 二 代 卵 的 密度 进行 了 插 
值 ,生成 的 各 年 空间 分 布 图 如 图 2。 
2.3 交叉 验证 

对 生成 的 分 布地 图 进行 交叉 验证 ,结果 如 表 2 
所 示 -。 根据 该 表 可 以 看 出 : 各 年 的 平均 预测 误差 和 
平均 标准 差 都 很 小 ,接近 于 0; 均 方 根 标准 差 也 接近 
于 1; 均 方 根 误 差 和 平均 克 立 格 标 准 差 均 相 对 较 小 。 
因此 各 年 的 模拟 分 布 图 基本 达到 了 比较 小 的 误差 ， 
得 到 了 较 好 的 模拟 效果 . 
2.4 ” 卵 量 的 了 时 间 序 列 分 析 

为 了 对 研究 区 1980 ~ 1997 年 18 年 的 发 生 情况 
有 一 个 整体 了 解 , 首 先 对 各 年 所 有 抽样 点 的 密度 值 


30 昆虫 学 报 Acta Entomologica Sinica 


表 1 各 年 的 插值 方法 和 相应 的 模型 参数 
Table 1 JInterpolation methods and parameters of the models each year 
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空间 变异 成 分 ”随机 程度 


相关 程 和 太 使 
插 信 方法 年份 。 模型 内 变量 Be Me 
Interpolation | : Sill Partial sil Nugget . 
a Year Model Covariable the spatial (EY CC EY Le) Proportion of degree 
me Ss dependence (km) " y spatial variation 9 To 
6 月 上 名 降雨 量 
Ordinary Cokriging ”1980 球 型 Spherical Gross rainfall in the 115.95 39984 16080 23 904 40.22 59.78 
first ten days of the June 
棉花 种 植 比例 
Ordinary Cokriging ”1981 ” 球 型 Spherical Proportion of the 283.52 12408.3 439.29 11 969 3.34 96.46 
cotton planting 
Ordinary Cokriging 1982 ”高 斯 Gaussian 高 程 levation 293.66 448 800 153 400 295 400 34.18 65.82 
5 月 降雨 强度 
Ordinary Cokriging ” 1983 ”高 斯 Gaussian Rainfall intensity 189.13 149623 60970 88 653 40.75 59.25 
in May 
Simple kriging 1984 ”指数 Exponential 363.66 63234 10433 52 801 16.50 83.50 
Ordinary kriging ”1985 ”指数 Exponential 473.67 32295.6 4552.6 27743 14.10 85.90 
6 月 降雨 总 量 
Ordinary Cokriging 1986 ”高 斯 Gaussian Gross rainfall 480.33 197 437 139337 181 900 7.88 92.12 
in June 
Ordinary kriging ”1987 ” 球 型 Spherical 317.59 11227.4 2546.8 83680.6 22.68 77.32 
Ordinary kriging ”1988 ” 球 型 Spherical 379.3 115706 49140 66 566 42.47 57.53 
Ordinary Cokriging ”1989 ”高 斯 Gaussian 高 程 有 Levation 185.12 146 852.4 7092.4 139 760 4.83 95.17 
Ordinary krigin 1990 ”指数 Exponential ， 343.74 32728.1 4782.1 27 946 14.61 85.39 
9 5 月 温 雨 系数 
二 237.16 435550 268 590 166 960 61.67 38 .33 
gr0sS Iainfa]] to average 
a temperature in May 
Ordinary kriging 1992 ”高 斯 Gaussian 487.08 8781 700 4 728 600 4 053 100 53.85 46.15 
Ordinary kriging 1993 ”高 斯 Gaussian 272.63 4222 910 114 510 4 108 400 yl 97.29 
Ordinary kriging 1994 ” 球 型 Spherical 417.71 90447 34635 C55 812 38.29 61.71 
Ordinary Cokriging 1995 ”高 斯 Gaussian 358.34 2515 100 1 205 300 1 309 800 47.92 52.08 
Simple kriging 1996 ” 球 型 Spherical 49.065 951960 416090 535 870 43.71 56.29 
5 月 下 旬 降 雨量 
Ordinary Cokriging ”1997 球 型 Spherical Gross rainfall in the last 154.91 411 154.3 74.348 411 080 0.02 99.98 


ten days of May 


取 平 均值 ,做 时 间 序 列 分 析 如 图 3。 结 果 表 明 18 年 
间 , 棉 铃 忠 卵 的 发 生 数 量 始终 处 于 波动 状态 , 其 中 
1980 ~ 1991 年 维持 在 低 水 平 的 流动 状态 , 其 发 生 的 
平均 密度 为 330 粒 / 百 株 , 变化 范围 为 122 ~ 706 镁 / 
百 株 ; 而 1992 ~ 1997 年 的 发 生 数 量 在 高 位 波动 , 其 
平均 密度 为 1 877 粒 / 百 株 ,变化 范围 为 264 ~ 4 023 
粒 / 百 株 。 其 中 1992 年 发 生 数 量 最 高 , 该 年 平均 值 
为 4 023 粒 / 百 株 , 远 远 高 于 其 它 年 份 , 其 次 为 1993 
年 ,发 生 密 度 为 2 820 粒 / 百 株 。 从 整体 的 发 生 情 况 
看 ,1992 年 发 生 最 重 , 其 次 为 1993 年 和 1995 年 
C2 178 粒 / 百 株 ), 1981 年 (145 粒 / 百 株 ) 和 1990 年 
C144 粒 / 百 株 ) 发 生 较 轻 , 1987 年 发 生 最 轻 (122 粒 / 
百 株 ),1980 年 代 的 发 生 数 量 显 著 低 于 1990 年 代 。 
2.5 空间 分 布 及 动态 

根据 空间 相关 分 析 的 结果 (各 年 卵 在 空间 均 存 
在 空间 相关 性 ,并 具有 一 定 的 相关 程 ) 以 及 各 年 模拟 
地 图 都 存在 聚集 斑 块 等 特点 ,可 以 确定 研究 区 二 代 


棉铃 虫 卵 在 空间 呈 聚 集 分 布 的 格局 , 而 各 年 卵 的 聚 
集 程度 依据 其 空间 相关 成 分 所 占 比 例 的 不 同 ( 表 1) 
而 不 同 。 如 1991 年 的 空间 变异 成 分 所 占 比 例 在 各 
年 中 最 高 (61.67%), 其 卵 的 聚集 程度 亦 最 强 , 方 回 
性 也 最 明显 。 

从 空间 模拟 地 图 还 可 以 确定 其 在 地 域 上 的 发 生 
趋势 , 即 从 东北 至 西南 , 发 生 程度 逐渐 趋 重 ,其 中 西 
部 的 旱 平 县 .曲阳 县 . 行 唐 县 ` 灵 寿 县 .平山 县 等 地 以 
及 南部 的 邯郸 市 、 磁 县 、 临 漳 县 、 魏 县 等 地 发 生 最 重 ， 
而 东北 的 遵 化 市 、 玉 田 县 .丰润 县 .唐山 市 、 丰 南 县 等 
地 最 轻 。 

对 各 个 发 生 等 级 所 占 的 面积 和 比例 进行 统计 
( 表 3), 可 以 看 出 ,各 个 发 生 等 级 不 仅 在 地 域 上 存在 
差异 ,而 且 随 时 间 发 生变 化 : 1) 轻 发 生 等 级 仅 发 生 
在 1980 和 1987 年 ,但 所 占 面积 很 小 ,分 别 仅 占 当年 
发 生 总 面积 的 0.02% 和 1.77 多 :2) 中 全 轻 的 所 占 面 
积 在 大 部 分 年 份 都 比较 低 ,只 有 1981 年 CS$4.65 狗 人) 
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图 2 河北 省 二 代 棉 铃 虫 卵 的 空间 分 布 及 动态 
Fig.2 Spatial distrbution and dynamics of the 2nd generation eggs of the Helicoverpa armigera(l Hiibner) in Hebei Province 


粒 / 百 株 Eggs/100 plants 
多 SE 
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年 份 Year 


图 3 河北 省 二 代 棉 铃 虫 卵 的 平均 数量 动态 (1980 ~ 1997) 
Fig.3 Average quantity dynamics of the 2nd generation eggs of the Helicoverpa armigera( Hiibner) in Hebei Province (1980 - 1997) 


1987 年 (86.08% ) 和 1990 年 (76.53%) 三 年 所 占 面 
积 比 例 较 大 , 均 超 过 50% ;3) 中 等 发 生 等 级 所 占 面 
积 比 例 在 50% 以 上 的 年 份 有 6 年 ,分 别 为 1980 年 
(78.8% ) .1985 年 (89.24%)、1986 年 (63.7% )、1988 
年 ( 64.41%)、1989 年 (75.98%) 和 1994 年 
(99.61% );4) 中 偏重 等 级 发 生 比 例 在 50% 以 上 的 年 
份 有 4 年 ,分 别 为 1982 年 (72.09%)、1984 年 


(100 和 % ) 1996 年 (96.06% ) 和 1997 年 (97.92% );5) 
大 发 生 等 级 仅 有 6 年 存在 ,其 中 发 生 面积 比例 在 
50 以 上 的 年 份 有 3 年 ,分 别 为 1992 年 (95.85% ) 、 
1993 年 (94.47% ) 和 1995 年 (80.41% ) ;6)1983 年 和 
1991 年 各 个 等 级 发 生 面积 比例 均 未 达到 50% ,但 这 
两 年 的 中 发 生 等 级 所 占 面 积 比例 均 很 大 。 
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插值 方法 Interpolation methods 年 份 Year 


Ordinary Cokriging 
Ordinary Cokriging 
Ordinary Cokriging 
Ordinary Cokriging 
Simple kriging 
Ordinary kriging 
Ordinary Cokriging 
Ordinary kriging 
Ordinary kriging 
Ordinary Cokriging 
Ordinary kriging 
Ordinary Cokriging 
Ordinary kriging 
Ordinary kripging 
Ordinary kriging 
Ordinary Cokriging 
Simple kriging 
Ordinary Cokriging 


<50( 轻 发 生 Lowest) 


1980 
1981 
1982 
1983 
1984 
1983 
1986 
1987 
1988 
1989 
1990 
1901 
len 
1993 
1994 
1993 
1996 
1997 


51 ~ 150( 中 偏 轻 Lower) 
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表 2 各 年 插值 地 图 的 交叉 验证 结果 


Mean 
-0.030 
0.007 
0.039 
0.060 
0.039 
0.011 
0.003 
— 0.051 
0.487 
— 0.030 
0.091 
— 0.064 
0.078 
— 0.043 
— 0.056 
0.048 
0.042 
0.605 


RMS 
188.7 
103 .3 
363 .9 
283 .3 
234 .1 
161.7 
413.1 
103 
271.8 
362.2 
Ll 
472.3 
2 192 
1 983 
289.7 
1 300 
5968 .3 
637 .9 


151 ~ 500C 中 发 生 Moderate) 501 ~ 1 500( 中 偏重 Higher) 


Table 2 The results of cross validation for interpolation maps each year 


Avg. Std Mean Std RMS Std 


194.2 
114.2 
388.3 
339.8 
243.2 
LE 
443.6 
Rs 
202.8 
393.4 
180.3 
324.1 
2.108 
2 107 
2023 
1 269 
969.1 
672.8 


— 0.006 
0.0003 

— 0.0034 
— 0.0048 
— 0.0003 
—0.0014 
— 0.0003 
— 0.0006 
0.0008 

— 0.0023 
—0.0009 
-0.0010 
-0.0024 
— 0.0003 
一 0.0033 
-0.0010 
0.0009 

0.00003 


0.997 
0.923 
0.968 
0.846 
L087 
0.914 
0.935 
1.032 
0.939 
0.926 
1.057 
0.922 
1303 
0.944 
1.097 
让 
0.998 
0.978 


RMS 
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— Avg. Std 
-5.5 
-9.2 


-24.6 
一 34.3 


8.9 


-16.2 
-28.35 


| 
| 


= 人 S12 


10.8 


29 
一 122 
1 :2 

Sl 
-0.6 


-14.9 


Mean: 平均 预测 误差 Mean prediction error; RMS: 均 方 根 误差 Root-mean-square predicition error; Avg. Std: 平均 克 立 格 标 准 差 Average kriging 
standard error; Mean Std: 平均 标准 差 Mean standardized prediction error; RMS Std: 均 方 根 标 准 差 Root-mean-square standardized prediction error; RMS - 
Avg. Std: 均 方 根 误差 与 平均 克 立 格 标准 差 之 间 的 差 Difference between Root-mean-square prediction error and Average kriging standard error. 


表 3 各 发 生 等 级 的 面积 和 比例 (1980 ~ 1997) 
Table 3 Areas and proportions of the outbreak grades (1980 - 1997) 


> 1 500C 


重 发 生 Highest) 


Year 面积 Area 比例 Proportion 面积 Area 比例 Proportion 加 积 Area 比例 Proportion 面积 Area 比例 Proportion 面积 Area 比例 Proportion 
(1Lx1l0 km)  (%) 


(1x 10 km2) 
1980 0.0017 


1981 
1982 
1983 
1984 
1985 
1986 
1987 0.1659 
1988 
1989 
1990 
1991 
1992 
1993 
1994 
1993 
1996 
1997 


(的 ) 
0.02 


1.77 


1.982 
3.1123 
0.3985 

1.393 


1.0069 
3.3959 
8.0533 
0.906 

0.9713 
7.1391 
0.6144 


0.0366 


21.19 
34.63 
6 .4 
17.03 


10.76 
36.30 
86.08 
9.68 
10.38 
76.33 
6.37 


0.39 


(1 x 10 km ) (%) 


7.3714 
4.2426 
2.0129 
4.5367 


8.3482 
3.9392 
1 .1359 
68.026 
7.108 
2.196 
4.2292 


9.3185 
0.5174 


0.1823 


78.80 
45.33 
21.32 
48.49 


89 .24 
63.70 
12.14 
84.41 
75.98 
23.47 
45.21 


99.61 
33 


1.73 


(1 x 10 km ) 


6.7437 
3.2234 
9.3531 


之 .4231 
1 .27358 


3.6691 
0.388 
0.3176 


1.3148 
8.9865 
9.1606 


(WW) 


72.09 
34.48 
100.00 


25.90 
13.64 


39.22 
4.15 
3.33 


14.03 
96.06 
97 .92 


(1x 10 ko2) (%) 
0.8424 9.00 
8.9671 95 .85 
8.8375 94.47 
7.5229 80.41 
0.3686 3.94 
0.0322 0.34 


基于 以 上 结果 ,车 根据 等 级 发 生 面 积 比 例 大 于 年 ;中 发 生年 包括 1980、1985、1986、1988、1989 和 


50% 作 为 衡量 各 年 发 生 程 度 的 指标 , 则 可 以 将 各 年 
的 发 生 大 致 划 分 为 轻 发 生年 , 如 1981、1987 和 1990 


1994 年 ; 大 发 生年 包括 1982、1984、1996、1997、1992、 
1993 和 1995 年 ;在 这 18 年 当中 , 1992 年 发 生 最 重 ， 
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其 次 为 1993 年 和 1995 年 ,1987 年 发 生 最 轻 。 
3 讨论 


害虫 一 般 都 滋生 于 一 定 地 域 的 寄主 作物 或 植被 
之 中 ,存在 于 某 一 特定 的 时 间 段 内 , 也 就 是 说 它 的 生 
态 学 变量 既 有 时 间 性 ,也 有 空间 性 。 比 如 旷 的 密度 
数值 为 $0 粒 / 百 株 ,其 生态 学 意义 就 是 某 一 时 刻 (如 
1980 年 6 月 18 日 上 午 10 点 )、 某 一 地 点 (具有 确定 
的 经 纬度 值 ), 每 一 百 株 棉花 所 具有 的 卵 粒 数 。 从 本 
文 的 空间 结构 分 析 来 看 , 尽管 结构 性 成 分 所 占 比 例 
不 是 很 大 ,但 这 个 结构 性 成 分 在 各 年 确实 都 存在 , 它 
是 区 域 化 变量 的 组 成 部 分 ,不 应 被 忽略 。 

尽管 空间 相关 性 是 客观 存在 的 ,但 各 年 的 相关 
距离 (相关 程 ) 却 存在 差异 ,空间 变异 成 分 所 占 比例 
普遍 较 低 。 在 相同 的 研究 区 域 出 现 不 同 的 相关 距离 
是 由 空间 异 质 性 在 各 年 的 差异 所 导致 的 ,这 种 空间 
异 质 性 的 差异 包括 寄主 作物 的 空间 分 布 .气象 条 件 
的 不 同等 等 ; 而 空间 变异 成 分 所 占 比例 偏 低 《 相 对 于 
随机 性 成 分 ?是 由 两 方面 的 因素 导致 的 : 一 是 各 年 
的 样 点 数 多 为 30 个 左右 ,参与 插值 的 样 点 数 偏 少 。 
一 般 而 言 , 参 与 插值 的 样 点 数 越 多 佑 值 越 准确 , 越 少 
误差 越 大 。 二 是 各 地 发 生 数量 差异 较 大 , 数据 方差 
较 大 , 随机 性 增 大 。 但 从 结果 仍 可 以 看 到 空间 结构 
成 分 在 各 年 仍 占 有 相当 的 比例 , 其 中 最 高 的 有 
61.67%(1991 年 )。 

在 空间 相关 性 分 析 的 基础 上 , 我 们 模拟 了 卵 在 
各 年 的 空间 分 布 ,这 种 模拟 不 是 简单 的 距离 加 权 , 而 
是 建立 在 对 整个 研究 区 所 有 抽样 点 的 空间 相关 分 析 
基础 上 的 。 不 仅 对 单一 变量 进行 简单 的 元 立 格 插 
值 ,也 应 用 了 协同 克 立 格 方法 , 它 综合 考虑 了 影响 旷 
量变 化 的 其 它 变 量 并 比较 了 各 种 方法 和 模型 的 基础 
上 对 卵 的 密度 进行 模拟 , 因而 最 大 限度 的 降低 了 模 
拟 的 误差 ,提高 了 模拟 的 精确 性 。 

就 本 项 研究 而 言 ,通过 对 18 年 的 卵 的 密度 数据 
的 分 布 模拟 , 我 们 可 以 从 时 空 两 个 角度 直观 地 分 析 
它 的 动态 变化 , 即 可 以 定性 地 判断 和 定量 地 分 析 并 
确定 卵 的 数量 在 不 同年 份 的 发 生 程度 及 变化 , 确定 
同一 年 份 卵 在 空间 是 如 何 分 布 的 , 重 发 生地 位 于 何 
处 。 同 时 间 序 列 相 比 ,空间 分 布 及 变化 的 研究 不 仅 
能 够 说 明 时 间 序 列 常 用 的 平均 的 数量 动态 , 还 提供 
了 空间 连续 的 数量 分 布 及 详细 的 地 理 环 境 信 息 。 但 
空间 分 布 的 模拟 并 不 只 是 了 解 一 下 害虫 的 动态 信 
县 ,更 重要 的 是 揭示 害虫 时 衬 动 态 的 内 在 机 制 ,并 为 


害 忠 的 宏观 管理 和 预警 莫 定 基础 。 因 此 ,还 需要 结 
合 地 理 和 和 气象 等 因子 对 害虫 时 空 分 布 及 变化 的 穆 因 
作 进一步 的 分 析 。 

就 棉铃 虫 而 言 ,影响 其 种 群 动态 的 因子 很 多 , 其 
中 包括 地 理 的 如 所 在 区 域 的 地 理 位 置 .地 势 ,气象 的 
如 温度 和 降雨 ,植被 覆盖 如 寄主 作物 的 种 类 及 其 种 
植 比例 ,人 为 因素 如 耕作 制度 等 , 这 些 因子 在 空间 的 
分 布 及 变化 将 直接 或 间接 地 影响 棉铃 忠 的 种 群 数 
量 ,或 者 说 棉铃 忠 种 群 数量 在 时 间 和 空间 上 的 变化 
是 这 些 因 子 综合 作用 的 结果 。 

许多 研究 表明 : 虫 源 基 数 、 食 料 和 气候 条 件 是 
影响 棉铃 虫 种 群 动态 的 三 个 主要 因素 ,在 虫 源 基 数 、 
食料 相对 稳定 的 条 件 下 ,气候 条 件 即 温度 和 降雨 成 
为 影响 其 发 生 的 关键 ( 华 芜 楠 , 1995; 华 医 楠 等 ， 
1996)。 株 铃 虫 在 1990 年 代 的 发 生 数 量 显 闭 高 于 
1980 年 代 , 其 主要 原因 就 是 气候 条 件 的 变化 ( 夏 获 
源 ,1994; 秦 淑 莲 等 ,2003)。 棉 铃 虫 适 生 于 半 干 旱地 
区 ,成 虫 繁殖 喜 高 温 高 湿 , 而 进入 1990 年 代 以 来 , 冬 
季 普 遍 偏 暖 ,翌年 春季 和 气温 回升 快 。 气 温 升 高 缩短 
了 各 代 虫 的 发 育 历 期 ,使 各 代 发 生 期 提前 ,暴发 机 率 
增 大 ( 谢 宝 瑜 等 , 1996)。 

研究 区 棉铃 虫 重 发 生地 位 于 西部 及 西南 部 分 地 
区 ,这些 地 区 的 主要 特点 是 多 山 、 植 被 丰富 、 小 麦 种 
植 比例 高 以 及 6 月 降雨 量 偏 少 。 多 样 性 的 寄主 可 以 
为 柳 铃 忠 一 代 幼 忠 提 供 丰 富 的 食料 , 同时 也 提供 了 
避 护 场所 和 转移 为 害 的 条 件 ( 谢 宝 瑜 等 ,1996), 而 6 
月 降雨 量 依 少 则 有 利于 二 代 棉 铃 虫 的 生长 和 繁殖 ， 
从 而 也 容易 导致 二 代 在 栅 田 的 发 生 数 量 侦 重 。 
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